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TECHNICAL NOTE 2439 
/ 

A THEORY OF CONDUCTIVITy OF COLD-WORKED COPPER 
, Rolf Landauer 


SUMMARY 

The Increase in the resistivity of copper under cold-"WDrking is 
calculated. The increase is eissumed to he caused hy dislocations sur- 
rounded hy a long-range electrostatic field that _ scatters the conduc- 
tion electrons . The amount of scattering is found hy the method of 
deformation potentials of Bardeen and Shockley. The scattering is 
present in addition to' the normal thermal scattering and is regarded 
as a peidiurhation in the Boltzmann equation. This pertiirhation is used 
to find the incremental resistance per dislocation. From this calcu- 
lated increment in resistance and the known increase of resistivity of 
heavily cold-worked copper, the number of dislocations in the cold- 
worked copper is fomd to he in agreement with the number estimated on 
the hasis of stored-energy measurements. 


INTRODUCTION 

In the vicinity of an edge-type dislocation, a metal is strained. 

A dilation of the lattice is associated with this strain; therefore the 
density of electrons varies in the vicinity of the dislocation and the 
■\7ldth of the filled portion of the conduction hand must also vary. In 
equilihrlum, however, the top of the filled portion of the conduction 
hand must he at the same level everjadiere; so the bottom of the conduc- 
tion hand must accommodate the variation. This variation scat o. 
electrons and is taken into accovint in the ensuing calciilation. 

The matrix elements for this scattering are calculated hy the 
method of deformation potentials of Bardeen, and Shockley (reference l). 
Once these matrix elements have been obtained, the rest of the calcu- 
lation foUotra the method of Meckenzie and Sondhelmer (reference 2), 
•vdiich treats the scattering due to dislocations as a small perturbation 
in the Boltzmann eqxiatlon. This treatment gives the change of resist- 
ance at high temperatures, idaere the increment is only a few percent 
of the normal thermal resistance. An earlier calculation hy Koehler 
(reference 3) deals with the resistance , of cold-worked copper at abso- 
lute zero. / 
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The physical mo'biva’bion. for this calculation and the results have 
"been published in reference 4, "whereas the detailed procedure is given 
in this report. 


AHALYSIS 

In order to discuss electronic motion in the "vicinity of an edge 
dislocation, a coordina"te system must he defined. The z-axis is taken 
"to he the dislocation axis and "the distance of a point from this axis 
is denoted hy r^. The angle of elevation above "the slip plane is 

measured hy 0 . The dislocation "will he taken "to he positive, so that 
there "will he an extra plane of a"boms above "the slip plane at 0 = «/2. 
The slip direction 0 =* 0 is considered parallel to "the x-axis. (Sym- 
bols aie listed in appendix A.) 

The density of ions varies near the dislocation. If iiq is the 
density -of ions in "bhe "unstrained metal and if An is "the increment in 
the number of ions per unit volume "then, as sho^ in appendix B, 

^ a (1-P-2P^) sin 0 

ri, 

where a is "the slip distance, and V is Poisson’s ratio. The "width 
of the filled portion of the conduction hand occupied hy n electrons 
"wi"th effective mass m* is 



Kot only n, hut also m* is a function of position in "the neighhorhood- 
of the dislocation. 

TOien "two different metals are hro"ught into contact, an electro- 
static field is set "up so "that "their Fermi levels are brought together. 
The hot"tom levels of "the "two conduction hands are "then at different 
energies, the difference being equal to the difference in the "wid"th of 
the conduction hands. In a dislocation, regions of the metal siahject 
to different strains are in contact "wi"th each other j. in the same "way as 
in the case of different metals it can then he expe^ed "that the Fermi 
levels are hrou^t "together in the vicinity of the dislocation, and that 
the variation in hand •wld"th gives rise "to a "variation in the energy of 
"the lo"wer hand edge. 

The method of deformation potentials (reference l) describes the 
most general me"thod hy "which electronic motion in deformed crys"tals may 
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"be treated. Let "be the wave fimction at the lower hand edge. Let 

5U(r) he thh deviation of -the lower hand edge from its normal position. 
Then the wave function for an electron with energy E is given hy 
A(r) ijfQCr) where A satisfies 

(3) 


(4) 

where is the effective mass in the -undistorted lattice. Since 

i^-Q^r) is taken to he the same for all wave functions in the hand, it 
need not he considered and only the perturhation in A(r) due to the 
dislocation must he evaluated. 


or 


ij- 2 Am 

+ 5U(r)| 

2m*(r) 


A(r) = EA(r) 


[2mo* mo* J 


r) = 0 


Now, for an electron near the top of the filled portion of the 
conduction hand, E-8U is the iriLdth of the conduction hand given hy 
equation (2); hence 



(E-6U) = 



Setting n =■ riQ + An gives, to the first order in An 

b3q ohiq \oit y o \p3t / 1 


(5) 


( 6 ) 


The second texm on the ri^t-hand side represents a perturhation due to 
the elastic distortion of the metal. Using equation (l) gives the per- 
turhation as 

2 r h^ _a^ (l-u-2P^) sin 9 

2 |_2mo*W/ J 2 rt i?) r^ 

or 

2 a (1-P-2P^) sin 9 

3 "^^0 2 % r^ 


( 8 ) 
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■viiiere Eb^O "widtli of the filled portion of the conduction hand 
in the unstrained metal. Therefore A(r) is to he , determined from the 
equation 


v2 + E_ . + E„ „ i liz!i:2£l sisi] A(?) 


The unpenrburhed equation is 


and its solution is AqC?) = ^ ■where 


(9) 



®B,0 


0 


and Ao(r) is the factor that modulates iJfo(r)j for an electron at the 

top of the conduction hand, in the undeformed i^etal. The prohahllity of 
transitions to states of the type 

A(?) = e^^*" 

must no'^f he calculated. The new state has the same energy as "the Ini - 
tial state 

6g|k- |g 

2mo 

httt represents the scattered electron traveling in a new direction. 

From this sca'ttering prohahility,. the increase in resistance cem he 
found hy the me'thod of reference 2. 

Note that it is not necessary to know m»(r), h\jt that it is only 
the -variation of n(r) that mat-tersj m*(r) need not he known. 

Equation (9) is the wave equation that a free particle of mass Po*^ 
incident energy charge -q oheys,. if scattered hy an elec- 

tros-fcatic potential 


^,0 
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V = ^ (l-l?-2P^) sin 6 

^ 2ng r^. 

This Is the potential of a line dipole located along the dislocation 
axis . 


CALCULATIOIT OF MATRIX ELMENTS 
The perturbing potential in equation (9) has the form 

V = p SiSJ. ( 10 ) 

I’t 


where 


2 _a^ (l-l?-2l^) 

3 2it 


( 1 -^) 


%,0 


( 11 ) 


!This potentleil is assumed to exist in a rectangular box of length L 
along the dislocation axis and of area d^ measured in a cross section 
perpendicular to the length L of the box. The unperturbed wave func- 
tions are. normalized for this box, are selected to satisfy periodic 
boundary conditions, and are of the form 


Ao(r) = - 7 = e' 


it*r 


( 12 ) 


where V = Ld and is the volume of the box. 
The matrix element of interest is 


1 

V 


:>i(k’-S)»r 


p siE-§.dY 
^t 


(13) 


■which vanishes \inless = kg’. Let k^, k|.’, and r.^ represeht the 

congponents of k, k’, and r, respectively, that are perpendicular to 
the dislocation axis. After integration along the z-axLs, e^^res-' 
Sion (13) becomes (if kg = k'g) 


f 
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i(k. »-k. ).,r 




sin 0 
"•t 


d\ 


(14) 


•where d\. is an element of surfane of ^ p^ane perpendicular to the dis- 
location axis. Let the magnitude of he denoted hy K and 

the angle it makes ■with the slip plane hy 
in question becomes 

iKr^ cos (S-Sq) 
e P 

Setting 0-0Q = ^ gives 



0Q. Then the matrix element 

* 

5%^ aK (15) 


^1' 


iKr+ cos a o 

e — sin (a+e^) dX 

X rt ' 0' 


(16) 


•where sin (o/+5q) can he escpanded as 


sin a cos 6q + cos a sin Gq 


(17) 


The first term of expression (17) is odd in a. Because the exponential 
is even in a, this contribution to the integral vanishes. The 
remaining integrand gives 


P 


P 

^0 iKr-t cos a 
d^ 


(18) 


The integration over r^ from r-jj = 0 to r^ = ot> can he performed hy 
assuming that k has a small positive imaginary part and then letting 
this imaginary part approach zero. This procedure gives for the matrix 
element 


.».» P sin 9n I 1 
M(k,k’ ) = — da 

d^ J 

2n p sin Gq 
iK d^ 


(19) 


( 20 ) 
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Let 0^‘ and 0^ te the angles that "Sj.’ and Icj. (final and 
initial wave vectors, respectively) make with the slip plane. Then 



If jlt'i - |tt|. 



Now consider this transition in kj. space, as shown in this fig- 
ure. The plane represents all the states i^lth a given value of k^. 

The curve represents a contour of constant energy. As is usual in cal- 
•culations of this type, the only transitions that are important are 
those that almost conserve energy. The electron is assumed in an ini- 
tial state specified by k^ . . The quantity needed is the transition 
probability to the. set of states in the element dS d|. The reasoning 
used here follows that given in reference 5. 

If is the energy of the final state and Ej^ is the energy 

of the ^^oitial state, the probability that the electron will be in 
state k’ after time t is 
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( 22 ) 


■5-ihere p = (Ej^* -I^)/li. How the number of states per unit area of this 

tTO-dimensional hj. space is — and the prohahillty of finding the 

d? 

electron in dS d| is therefore — ^ total proha- 

4jt 

hility for finding the electron in dS is 



(23) 


The matrix element M(^,k’) can he. assumed constant over the range of 
d S. The only important contrihutions to the Integrand come from small 
values of i for which 


P 


1 2E t 


(24) 


where i is taken to he zero on the contour of constant energy shown 
in the preceding figure. The quantity actually needed is the 
prohahillty of transition per unit time. This differentiation 
then gives 


P(k,k») dS = 



(25) 


Now this prohahillty of transition is superposed on the transition 
prohahillty arising from the thermal vihrations of the lattice. In the 
normal metal, in the annealed state, and in the presence of an electric 
field, there is a deviation g(^ from the Fermi distribution. This 
deviation is limited hy the thermal scattering, the amomt of the devia- 
tion determining the amoxint of the current. Cold working changes the 
conductivity of nonporoxis and isotropically conducting metals hy only a 
few percent. In the presence of the field, the dislocation, and the 
thermal scattering, the deviation from the Fermi distrihutlon is of the 
form g(k) - gj^(^ , where gj^ is a small perturbation and g(k) is 
the deviation that exists without the dislocation. Let foC^k) stand 

for the probability that the state associated with wave vector k is 
occiq)ied according to the Fermi distribution: 


2224 
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fo® = ^ 

1 + e 

t 


1 

[^(k) -^] /M? 


(26) 


Let f(k) stand for the actual distribution, f(k) = fo(k) + g(k) - g^^Ck) . 

The rate at which f (k) changes because of an electric field in the 
X direction is 


^ ^ ^ 

dt “ dk k -H 


I'There F is the field, an^ Tdiere the stirf aces of constant energy are 
assumed to be spheres ih k space. Furthermore, the conduction process 
in the annealed metsil will be assumed describable by a relaxation time 
time T (references 2 and 5) so that 


g(k) 


dk k -tt 3k^ kt •& 


(28) 


In equilibrium, 3f/^ = 0. The rate of change of f due to the 
field as given by equation (27), must therefore be balanced by the scat- 
tering. The rate of scattering due to the thermal vibrations is given 
ty 


df (f-fo) g ^ 81 

3t " ’ T " " T T 


and the rate of scattering due to the dislocation is 

■’) dS' 


or 


=J^ ^(^) - f(k^ P(k,k» 


[g(S>) - g(k)J P(k,k‘) dS* 

- P(k,k’ ) dS* 


(29) 


(30) 


( 31 ) 
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la each, inbegral occurring ia eq.xjation (31) fhe integration is over the 
circle defined ty Ki=KI , -which lies in the plane kg.’ = kg.. 

The first integral -vanlBhes hecause last 

Integral is neglected. Both gi(k) and p(?,k*) are proportional to 
the perturbation represented hy the dislocation. The last integral is 
therefore of second, order in the perturbation. The total rate of 
change is 


dk 


+ [g(?) - g(k)]p(k,k«) dS’ = 0 


(32) 


According to eq.uation (28) the first two terms cancel. The remaining 
equation gives 


-V 

= - T [g(k’) -g(k)Jp(k,k’) dS» 


Si = 

•where g is given by eq-uation (28), 

Substituting from equation (28), the solution for g]^ becomes 


(33) 


Si(t) 




^0 ’IF I 

U 


(cos - COB 0-t,*) P(k’,k) dS* 


(34) 


The -mlue of P(k’,k) dS’ is given by equation (25), Substituting 
this -value in equation (34) gives 


g.(k) = t2 ^ ^ ^ pL i 

Bkfc H B fl2 3e 

^ 'Sk^vj 


Now dS ' = kij dS.{. ' . Therefore 


_ ^0 1 SF 27f32 1 


(cos - cos 9^ 


’) ^ / dS’ 


(35) 


(cos 0t " 0.^’) 


cos 


sin 




d-df 


(36) 


P*rom this change in electron distribution, the change in the current 
density can be determined: 
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= 


2 Sl(k) 
-- 3 - 


2 rt" 


The electron velocity is given hy 

1 
fi 


Vt = — -F^ cos 

3. «c T» 


and furthermoi^ 


(37) 


(38) 


dT^= sin y d^ d 0 ^ dk 


(39) 


where 7 is the angle that the k vector makes with the dislocation 
axis . Hence 


M 



gl(k) q SE 
^3 


cos 0 t ^2 sin 7 d 7 d 0 ^ dk 


(40) 


In order to evalmte AJ, the expression given for gp(k) in eq.ua- 
tlon (36) must he inserted in equation (40), This suhstitution gives 


-\ 

Aj £ F a 


5fr 


dk d/ d 9 i d0+'k^i eln 7 00s ft 
^ kt ^ 


(00s &b - COS ftt’) coa2^-ii-S-^ 

‘ 


3in2 


Hov? kj. = k sin 7 , Furthermore, 

-<r — sin 7 

^ dk k dk ' 


(41) 


(42) 


Hence 




3 Q ^ o fl2 dFo 

AJ=— g^T^-kgF dkdy d0^ dft^ 'k sin 7 cos 0 ^ 


(cos 0-t ®t') cos'll g j 

W) 


sin?\ 


( 43 ) 



12 


MCA TU 2439 


Only near tihe top of the filled portion of the hand is dfo/dk appreciably 
different from zero; hence it can he considered a negative Dirac S 
function because all other terms in the integrand vary slowly with k. 

The integration over Xt frcm y - 0 to y = it, and the integration 
over k, from k = 0 to k = or^ give 




The particular value of k associated -with the Fermi level is denoted 
by In appendix C the value of the integral is shown to be 


Therefore 





(45) 


and for the change in conductivity (absolute value is given) 
Kow if there are U dislocations in the area d2 

IajI sz 4: 3^ kf -E. 

= 4 3^. t2 r 2 kp 


(46) 


(47) 


where ^ is the number of dislocations per sq.uar>e centimeter. 


Substituting the valrie of p given in egLuation (11) gives 



2 

2 


l~U-2u^ 

_l-u2 ^ 




(48) 


The percentage change in conductivity is giveh by |ac|/o. This is 
also the percentage change’ in resistivity, if the change is small,. The 
value of a is ng.^T/ino* where m^* is the effective mass in the 
ui^erturbed lattice, Mnce 


lA2l = JmI = i J_ Er kp mn* 


(^ 9 ) 
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Eg.mtlon (50) gives the change in resistivity for ouirrent flow in 
the slip direction. 

If the current flow is along the y-azis, that is, 0 = «/2, a 
slightly different computation applies. In eq,uations (27), (28), aM 
( 32 ), Is^ must he replaced hy This substitution gives, in^ead of 

eq.uation (34), 

p 

g 3 ^(k)'tt t 2 ^ (sin 0;j. - sin 0^*) P(k*,k) dS‘ 

ard an equivalent raplaoemenfc in all subsequent equations. Equation (38) 
must be replaced by 

Instead of equation (44), the following relation is then obtained: 

(sin 0^ - sin 0-fc’) cos^ ^ 




-±i!T^Fk.r 

d2 


d0^ d0-jj* slh 0^ 


sin^ 


m 


and in appendix 0 the value of the Integral is shown to be which 

gives 


a 


Ap 4: e.^ fl-v-2v^\ T ^ 2 ... V ^ 


In the z-dlrection, along the dislocation axle, ,the conductivity 
remains unchanged by the dislocation. 


EESDITS 

In a material in which dislocations occur with equal probability 
at all orientations, the increase in resistance is the average increase 
for the X, y, eind z directions (reference 2). Hence for an 
Isotropically cold-worked metal 


Ap _ 16 a^/l-u-2u^V T TP 2 




If the values given for copper in refe3?ence 4 are substituted, 

Ae = 4.0x10-1^ ^ 

P 
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Setting tills result eq.ua! to the 2-percent change usually ohserved in 
copper gives 

J/ = SXlO^/cm^ 

It is estl.mted (reference 6) that heavily cold-'Horked. copper has 

6x10^ dislocations per square centinieter. This estimate -was based 
on stored-energy measurements. The agreement is excellent^ indeed 
considerably better than the assimiptione involved warrant. 

The possibility that some of the stored energy may be in the form 
of screw dislocations has been neglected. In a simple cubic crystal, 
these dislocations are not accompanied by volume dilations and probably 
cause less scattering than edge-tjrpe dislocations. Furthermore, it 
should be talien into account that dislocations in a face-centered cubic 
metal, such as copper, occur in the form of half -dislocations. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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APPENDIX A 


SYMBOLS 

0?he foUowijQg symbols axe used in this report: 

Aoff) solution of effective mss equation in undeformed crystal 

jAjj.‘|2 protahility that transition from to ^ has occurred 

a slip distance 

d -width of hoz (along z-axie) used in quantizing electron 

states 


E energy of electron in stationary state in perturbed 

periodic potential 

width of filled portion of corduction hand, function of 
position in vicinity of dislocation 


%,0 


, width of filled porbion of conduction hand in unstrained 
metal 




energies of final and initial s-fcate, respectively 
F electric field producing current flow 

fQ("k5 prohahility that state ^ is occupied according to Fermi 

distrihution 



prohahility that state k is occupied in the presence of 
electric field and dislocation 


s(^S 

g(*^ - gi(S 


deviation from Fermi distribution in noiml undlstoited 
lattice, in the presence of an electric field 

deviation from Fermi distribution in the presence of dis- 
location and eleobrlo field 


h Planck’s constant 

h/2rt 

AS decrement in current due to dislocation 

t 


wave vector of electTOn in initial s-ta-be 
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k» 






V 

m*- 


U 


mTB Vector of electron In final state 

magnitude of k at top of filled portion of conduction 
tand 

components of respective vectors that are perpendicular 
to dislocation axis 

components of respective vectors along dislocation eixis 

length of hox (along z-axis) used in quantizing electron 
states 

effective electronic mass in unstrained lattice 

effective electronic mass, function of position in vicinity 
of dislocation 

matrix element for transition from state k to state k* 
due to perturhatlon hy dislocation 

O 

number of dislocations crossing area d 

number of dislocations per square centimeter 


n 


iIq + An 


n 


number of ions per unit volume in unstrained metal 


An 


change in number of ions per imit volume due to strain, 
a function of position in vicinity of dislocation 


P(^,k?) 


prd^hlllty per unit time that a transition from k to 
k' occurs 


p 

q magnitude of electronic change 

? position vector defining location vith respect to an origin 

on dislocation axis 


as 

0 ? 


component of that is perpendio^^lar. to dislocation axis 

varlahle of integration in fcj. space, dS is along a 
contoxir of constant energy 

absolute temperature 


t 


time 
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u 

T 

7‘ 

x,y ,2 


a 

P 

r 

C 

e 


®0 


K 

X 

V 

ae 


p 

Ap 

a 

Aa 


deviation of lower "barBi edge from level at which it is 
when far away from dislocation 

1/2(0^ + 0t') 

l/2(0t - 0t’) 

electrostatic potential whose scattering effect is 
eguivalent to dislocation 

perturbing potential due to dislocation in wave equation 

expectation value of electron velocity oon^onente 

rectangular coordinates with z-asls along dislocation and 
z-axis in slip direction 

0-00 

2 a g-p-2i>^) 

3 2rt ^,0 

angle that k vector makes with dislocation axis 
Fermi level 

angle of elevation above slip plane, measured from dls- , 
location axle 

angle kj.-k|j* makes with slip plane 

angles that Sj.' and respectively, make with the 
x-axis 

magnitude of ’ 

rectangle of area d^ perpendicular to dislocation axis 
Poisson’s ratio 

variable of Integration in ^ space; d| is perpendicular 
to a oontorir of constat energy 

resistivity of nomal meteil 

' change in resistivity caused by severe cold working 

conductivity of normal metal 

change in conductivity caused by severe cold working 
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+0 


relasation time for oondtiotioix prooees 
volxnme element in ^ space 

vave ftmotion for electron at 'bottom of conduction "band 
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APPEHDIZ B 


C5HAUSE OF ION DEBSHY NEAE A D3SLOCATION 

Tlie stress distribution surrounding an edge-type dislocation has 
been given in reference 7 as 

pn = - D(sin 30 + 3 sin 0)/2r 

Pg2 = D(sln ,30 - sin 0)/2r 

Pl2- = D(oos 30 + cos 0)/2r 

P33 ” “ D(4u sin 0)/2r (Bl) 

. Pl3 = ®23 “ 0 

vhere D = Ga/2«(l-u), G being the shear modulus, and u is Polsson*s 
ratio. Now the strain tensor S, , bas diagonal berms given by 



where p is the p3?eesure, or the trace of the stress tensor, and e is 

the unit tensor. The dilation is C and is therefore given by 

i 



(B3) 


Substituting from eq.uatlon (Bl) in eg.uation (B3) gives 



a (l-u-2u^) sin 0 
2Jt (1.1^) r 


(B4) 


The change in ionic density An/n is the jasgatlve of equation (B4) and 
is given by 


An a (1»^?-2 p^) sin 0 
a “ at (i_p2) r 
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APPEKDIS C 


E7ALI]Aa?I01J OF HJESERAIS • 


Ttis ircbegral 


It 


(cos 0^ - cos 0^^) oob2 


d.0. d0. * ops 0. 

■c "C "C 




.v)”^ vJ 

is to "be eTOlxtated. Let 
and. 


sin2 




(Cl) 


I (0t + ®t*) = 


I (0t - ®t*) ■= ^ 


(C2) 


(C3) 


Tlie Integral given tlien becomes 
'rt Pusst-v 


du dvjcos^ (u+v) - cos (wv) cos (u-vj-~~~ 


■=0vj^'’’“rt 


sin V 


'^VrrO ^ 


+ 2 


vjT=-lt^ 


U=Jt+V 


du dv ^os2 (u+v) - cos (u+v) cos (u-v^ — 


U=-Jt-V 


sin V 


(C4) 


Tile factor 2 .represents the Jacotian of tlie transfoimation, Tlie range 
of integration is the sq.uaro shown in the following figure: 
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Tlie first integral ranges over the dashed triangle | the second integral 
over the I’emaining portion of the square , 0?he terms in hrackets can 
also he written 

^2 sin^ u sin^ v - 2 cos u cos v sin u sin (C5) 

The first term of expression (C5) is even in v; the seconi is odd; 

oos^ u/sin^ V is even in v, The3refore the first term of expres- 
sion (05) produces contributions of like sign and equal magnitude in 
the two Integrals of eqmtion (C4), The second tem of expression (C5) 
produces contrihutions that cancel. Therefore expression (C4) reduces to 


"V V=rt p 

' 

J’STirO JU=V-Jt 


sln2 u cos^ u du dv 


(06) ■ 


The integral over u is a standard integral^ listed in tables of 
integration. This integration leaves 


X 

I -g sin 4v + 2rt ' - Z'n i 
. v=0 


= -g’ sin 4v + 27tv - v2 


2^2 _ jf2 


which is the value of (01) 
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The integral 




rK 


-Jt vJ 


00s2^f5l) 

t\ IL: ' 


-jt 


d0^ 40t* ®t(®^ ®t - ®t’) /V ' a 


is also to he evaluated. The transfoimatlon given hy equations (C2) 
and (C3) must he iised again and gives 

n. 


It p u=lt-v 

v=0^ u=v-if 
0 /^rurot+v 


du. dv Isin^ (u.+v) - sin (u+v) sin (u-v^ 

*- J sin^ V 


/v=-rtO 


du dv |sin^ (u+v) - sin (u+v) sin (u-vj 


U=-It-V 


oos^ u 
sin^ V 


The terms in the hrackets can he written 


^2 oos2 u sin^ v + 2 sin u cos v cos u sin v^ 


(09) 


(CIO) 


As in expression (01), only the term even in v contributes, leaving 

tl=It-V 


8 


■'it p tl=It-'V 

vJv=0 0u®=T-It 


cos^ u du dv . 


(Oil) 


After integration over u, the integral is 

It 


1 >- 


6v - 4 sin^ ^ ~ ^ 


;v=0 
= 6itv - 3v^ 


v), 


= Sit^ - 3it^ = 3rt^ 


( 012 ) 


■which is the value of expression (08). 
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